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ABSTRACT 


Scale  effects  in  cavitating  flow  are  considered  for  the  so-called 
limited  cavitation  flow  regime.  The  effects  on  cavitation  scaling  of  nu- 
clei and  dissolved  air  in  ordinary  water  are  considered.  Previous  work 
by  the  author  is  summarised  and  a theoretical  study  is  made  to  gain  in- 
sight into  the  relationships  that  mutt  hold  between  th»  parameters  which 
affect  the  inception  of  cavitation.  A simplified  theory  gives  only  rough 
qualitative  agreemtm  with  experime" 

PART  I 

INITIAL  CONSIDERATIONS 

Introduction 

1* 

The  piesert  wo;’.  - .--quo?  to  the  author's  'eliminary  report 

on  the  effects  of  body  size  and  free  ./.ream  velocity  ••oon  cavitation  on  an 
immersed  body.  The  nomenclature  of  this  report  ’v< ) l follow  that  used 
in  the  preliminary  pape  • . 

It  is  useful  to  divide  the  flow  of  a li^  arouud  . Mid  body  into 

three  regime s^which  are  called  nol'icav/iii?.^  *"• l ca  vital*  cn, 

and  full  cavity  flow,  respects  ’•/.  Th-’  eo&l:  •:*;/  '*■  »«  '/hi.-., control  the 
behavior  of  flow3  in  the  uon caviv . -7  egime  ' » t-.  -sft  JU'*11  known  for 

some  time.  For  full  cavity  fl  a,  Rclrh*  *•>...  s:^7  ~ *hat  the  essen- 
tial quantity  fer  determining  the  flow  geometry  is  the  ' vacation  number, 
K = (pQ-  FyJ/'j  pvj^t  where  p^  and  Vq  are  the  free  stream  values  of  the 
static  pressure  and  velocity,  respecti  c,l\ . / ir  'a-  liquid  density. 

Reichardt  takes  for  pv  sunr.  cf  r'J  -a  a pressure*  *rtun  the  cavity. 

In  the  present  work,  pv  will  denote  only  the  liquid  vapci  pressure. 
However,  for  the  limited  cavitation  flow  regime  there  is  a definite  lack 
of  knowledge  concerning  the  effects  of  body  size,  flow  velocity,  and  dis- 
solved air  content  upon  the  development  of  cavitatiuu.  -n  fact,  it  is  now 
customary  to  employ  the  cavitation  number  K as  the  orly  significant 
parameter  for  describing  all  inviscid  cavitating  flews,  it  is  shown  in 
this  paper  by  both  theoretical  and  experimental  means  that  the  use  of 
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only  the  cavitation  number  for  limited  cavitation  is  an  unjustified  simplifi- 
cation. 

The  object  of  this  study  is  to  determine  the  behavior  of  limited  cavi- 
tation on  a given  body  shape  placed  in  a rectilinear  flow  of  constant  free 
stream  velocity.  Such  a flow  is  approximated  in  the  teat  section  of  a 
water  tunnel  where  the  free  stream  velocity,  free  stream  static  pressure, 
and  the  amount  of  air  dissolved  in  the  water  can  be  controlled  as  required. 
If  one  considers  only  this  elementary  flow  configuration,  very  simple  ex* 
perirr^nts*  can  I'C  performed.  For  such  experiments,  the  spatial  pres- 
sure distributions  on  the  body  surface  will  be  known.  From  these  known 
pressure  distributions  and  the  free  stream  velocity.-  pressure -tints 
relationships  can  be  calculated  for  a particle  moving  along  the  body  sur- 
face with  the  liquid.  From  such  pressure -time  functions,  one  can  study 
analytically  the  behavior  of  incipient  cavitation  with  changing  body  scales. 

In  this  report  such  calculations  are  made  and  compared  with  experiment. 

Possible  Origin  of  Cavitation 

For  ordinary  untreated  water  it  will  be  assumed  that  there  are  nuclei 
containing  air  or  water  vapor,  or  both,  which  are  stabilized  on  small, 
solid  particles  in  the  liquid.  It  is  held  that  boiling  or  cavitation  must  be 
initiated  from  such  nuclei,  since  their  absence  would  mean  that  very  large 
surface  tension  forces  must  be  overcome  before  cavitation  can  start.  The 
slight  tension  under  which  cavitation  normally  occurs  lends  plausibility  to 
arguments  favoring  the  existence  of  the  nuclei.  If  one  supposes  that  cavi- 
tation originates  from  such  small  nuclei,  then  it  must  take  them  an  ap- 
preciable time  to  grow  to  a macroscopic  size.  It  is  clear  that  this  growth 
time  must  be  dependent  upon  the  pressure  to  which  the  nuclei  are  subjected. 
If  for  a constant  free  stream  pressure  and  velocity  only  the  scale  of  the 
immersed  body  is  changed,  the  liquid  flow  may  show  corresponding  changes 
in  the  cavitation  due  to  changes  in  the  time  available  for  bubble  growth. 
Under  steady,  or  almost  steady,,  conditions,  ordinary  untreated  water  has 
a definite  boiling  point  and  has  no  appreciable  tensile  strength.  Such 
water  will,  however,  withstand  tensions  if  it  is  subjected  to  transient  low 
pressures  of  short  duration.  For  water  flowing  around  a solid  body,  such 
short  duration  transient  tensions  can  be  produced  on  the  water  as  it  flows 


*The  experiments  performed  for  this  study  are  described  in  Ref.  I. 
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by  the  body  surface.  It  is  possible  that  the  differences  in  the  transient  ten- 
sions produced  by  changes  in  body  length  and  free  stream  velocity  can  ac- 
count for  corresponding  differences  in  the  state  of  cavitation.  It  is  the 
study  of  nucleus  response  under  these  transient  pressure  conditions  that 
is  of  primary  concern  here. 

Summary  of  Previous  Results 

In  the  preliminary  report,  the  assumption  that  cavitation  starts  from 
small  nuclei  in  the  liquid  gi^es  conditions  for  dynamically  similar  bubble 
growths,  it  was  found  that  no  useful  seeing  laws  for  limi*td  cavitation  can 
be  obtained  from  such  similarity  arguments.  In  fact,  the  similarity  calcu- 
lations indicated  that  the  cavitation  number,  K,  is  not  the  significant  param- 
eter for  defining  the  state  of  the  flow  in  the  limited  cavitation  flow  regime, 
and  that  for  flw..s  of  the  same  liquid  around  geometrically  similar  bodies 
dynamically  similar  growths  of  individual  bubbles  cannot  be  expected  to 
obtain. 

The  effects  of  body  size  on  cavitation  were  then  investigated  experi- 
mentally for  two  families  of  simple  shapes  in  a steady  rectilinear  flow. 

The  first  experiments  were  performed  on  a family  of  12%  thick  symmetW 

cal  Joukowski  hydrofoils  of  constant  chord  which  spanned  the  High  Speed 

3 

Water  Tunnel  te«-‘  section  Next,  R.  W.  Kermccn  supplied  some  experi- 
mental results  for  a aerx^..  of  rignt  ci’v.  Jiar  cylinders  with  hemispherical 
noses  in  "\xially  syi.. metric  flow. 

The  important  results  oi  the  experiments  may  be  summarized  as 
follows:  First,  it  was  observed  in  connection  with  the  high-speed  motion 
picture  studies  that  different  cavitation  numbers  were  required  to  obtain 
similar  cavitation  development  on  bodies  of  different  size.  This  behavior 
was  ther.  confirmed  at  one  free  stream  velocity  for  a range  of  cavitaticr. 
numbers.  Second,  it  was  found  that  the  cavitation  number  for  the  in- 
ception of  cavitation  increased  with  free  stream  velocity  and,  in  certain 
instances,  for  constant  velocity,  the  incipient  cavitation  number  was  found 
to  increase  with  the  size  of  the  body. 
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PART  II  aoU 

THEORETICAL  INVESTIGATIONS 

DUi  £ 

Introduction  to  Theory 

The  experimental  results  summarized  above  show  that  the  cavitation 
number  for  incipient  cavitation  exhibits  systematic  changes  with  variations 
in  the  free -stream  velocity,  and  that  for  some  conditions  the  incipient 
cavitation  number  at  constant  free -stream  velocity  changes  with  the  scale 
of  the  body.  It  is  the  purpose  of  this  xe-Mp--  to  investigate  the  conditions 
for  incipient  cavitation  by  analytical  means  so  that  more  precise  ideas  can 
be  obtained  about  the  relationships  between  the  parameters  which  influence 
the  behavior  of  the  cavitation. 

The  question  of  how  a bubble  grows  from  its  original  small  nucleus 
to  a macroscopic  size  is  a dynamic  problem,  since  the  observed  varia- 
tions in  the  inception  of  cavitation  are  to  be  ascribed  to  differences  in  the 
transient  tensions  on  the  water  produced  by  changes  in  body  length  and  free 
stream  velocity.  The  equation  of  motion  for  a spherical  bubble  of  radius 
R(i)  in  an  unlimited  incompressible  inviscid  liquid  is 

rr  + |r2  = [p(R)  - P(t)]/r  , (l) 

where  p is  the  liquid  density,  p(R)  is  the  liquid  pressure  at  the  bubble 

wall,  and  P(t)  is  the  external  pressure  field  in  the  liquid  far  from  the 

bubble.  The  superscribed  dots,  as  usual,  stand  for  differentiations  with 

respect  to  time.  Let  tr  be  the  surface  tension  of  water,  pv  the  vapor 

pressure  corresponding  to  the  temperature  of  the  liquio,  and  pa  the 

initial  air  pressure  within  the  bubble  before  any  growth  from  its  initial 

size  R has  occurred.  If  '.he  bubble  is  assumed  to  expand  under  iso- 
o 

thermal  conditions,  then 

2<r  /Kn\ 

P(R)  ■ PV  - --  * *a(— ) • (2) 

P(t)  will  be  determined  by  the  flow  over  the  solid  body. 

Basic  Assumptions  and  Definitions 

Before  proceeding  with  the  detailed  formulation  of  the  theory,  a 
statement  of  the  basic  asrum.'tior.i  underlying  the  whole  analysis  is  re- 
quired. 
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First,  vi*3Cov*  effects  are  entirely  neglected  in  the  present  work. 

For  example,  it  is  known  that  the  minimum  pressures  in  the  liquid  occur 
on  the  immersed  body.  In  fact,  the  pressures  decrease  slightly  as  one 
proceeds  through  the  boundary  layer  toward  the  body.*  In  view  of  the  very 
imall  huCleus  size  it  would  seem  that  at  least  in  the  initial  stages  of  bub- 
ble growth  the  phenomena  might  be  restricted  entirely  to  the  boundary 
layer.  The  consideration  cf  boundary  layer  effects  requires  further  experi- 
mental information  which  is  now  being  accumulated.  The  present  discus- 
sion will  not  take  them  into  account. 

the  ^erection  between  the  flow  around  t’..  ~ ' nc'.ng  buhbie 

and  the  flow  of  the  liquid  around  the  immersed  boay  is  not  considered.9 
Thus,  the  motion  of  the  expanding  bubole  will  be  treated  as  though  the 
liquid  is  infinite  in  all  directions  and  the  velocity  and  pressure  relation- 
ships on  the  model  will  ignore  the  presence  of  cavitation.  The  effect  of 
the  flow  around  the  body  is  related  to  the  bubble  growth  through  the  P(t) 
term  in  the  equation  of  motion  for  the  bubble. 

Third,  the  bubble  is  assumed  to  move  with  t*e  fluid.  For  the  very 
small  bubbles  considered  here,  buoyant  forces  ar:  small  and  the  viscous 
drag  will  be  high  so  that  any  relative  motion  between  the  bubble  and  the 
water  will  be  very  small. 

Fourth,  the  bubbles  are  assumed  to  be  spherical. 

In  tins  stuav  jf  U that  cavitation  is  initiated  from  s ’.I 

fr 

nuclei  which  contain  air  or  water  vapor,  or  both,  stabilized  on  small, 
sol  d particles  in  the  liquid.  In  the  noncavitating  flow  regime  it  is  as- 
sumed that  the  nuclei  do  not  have  an  opportunity  to  grow  into  bubbles  of 
macroscopic  or  visible  size.  The  difference,  then,  between  the  no”- 
cavitating  and  cavitating  flow  regimes  is  that  in  the  latter,  the  nuclei  are 
exposed  to  a pressure  environment  favorable  to  bubble  growth  for  a period 
sufficient  to  allow  for  the  appearance  of  macroscopic  bubbles.  The  ap  - 
pearance  of  macroscopic  bubble:  depends  upon  the  response  of  the  nuclei 
to  the  transient  low  pressure  created  by  the  flow  of  water  f.round  the  body. 
This  bubble  growth  problem  will  be  determined  by  two  conditions.  First, 
the  initial  conditions  R(0)  = 0 and  R(0)  = Rq  will  be  prescribed. 

*It  seems  probable  that  the  nuclei  in  a liquid  have  a range  of  sites, 
but  it  will  be  assumed  here  that  all  nuclei  are  of  equal  effective  radius  R . 
Their  size  will  be  estimated  in  a subse  quent  section. 
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Second,  the  pressure -time  relationship  or  "forcing  function"  which  acts 
on  the  nuclei  and  promotes  their  growth  is  prescribed  if  one  knows  the 
pressure  distribution  around  the  immersed  body,  the  body  size,  the 
free  stream  velocity,  and  the  free  stream  static  pressure.  This  simple 
determination  of  the  forcing  function  is  a consequence  of  the  fundamental 
assumptions  listed  above. 

For  a given  flow  situation,  (he  time  available  for  bubble  growth  is 
fixed  by  the  size  of  the  solid  body  and  the  free  stream  velocity.  Incipient 
cavitation  •**■*»•  ce  said  to  ex: - 1 if  a nucleus  grow."  from  its  initial  radius 
Rq  to  a radius  of  one  millimeter  during  the  time  it  is  exposed  tc  the  low 
pressure  which  favors  growth.  This  value  for  the  final  bubble  radius 
was  selected  because  a bubble  of  one  millimeter  is  visible  to  the  unaided 
eye  and  this  arbitrary  bubble  size  is  of  the  same  order  of  magnitude  as 
the  bubbles  of  visible  incipient  cavitation  on  bodies  ir.  the  water  tunnel. 

Using  this  experimentally  derived  estimate,  one  can  now  ask:  Given  a 

fixed  time  for  growth,  determined  by  the  free  stream  velocity  and  body 
size,  what  will  be  the  free  stream  static  pressure  pQ  (and  hence  K)  at 
which  the  final  bubble  radius  will  be  one  millimeter? 

Stability  of  Gas  Nuclei 

For  practically  all  cases  of  technical  interest,  the  water  will  con- 
tain dissolved  air.  Accordingly,  the  nuclei  will  contain  air  ae  well  as 

- - 

water  vapor.  However,  if  surface  tthsion  forces  act  on  these  small 
bubbles,  the  air  will  be  driven  from  the  nuclei  into  solution  with  the  sur- 
rounding water.  ^ On  the  other  hand,  a stable  nucleus  can  exist  if  air  or 
vapor  bubbles  a*e  attached  to  small  solid  particles  in  the  liquid.  In  this 
cas?  the  effective  surface  tension  must  be  zero.  Therefore,  if  the  effective 
surface  tension  is  initially  zero,  the  nuclei  can  exist  indefinitely,  and  by 
Henry's  Law  the  initial  air  pressure  in  the  gas  pocket  is  proportional  to 
the  concentration  of  air  dissolved  in  the  water.  As  the  bubble  grows  from 
its  initial  effective  radius  R0,  the  surface  tension  will  increase  to  the  full 
value  2<r/ R.  The  required  behavior  of  the  surface  tension  may  be  approxi- 
mated by  p = 2S(R,  <r)/R  where  the  "surface  tension  law",  S(R,c  ),  is 
<r 

characterized  by  S(R0,<r)  = 0 and  SiR^.o1)  = <r  . Let  R^  be  that  bubble 
radius  at  which  the  surface  tension  law  first  achieves,  its  full  value  and 
put  R,  = nRQ.  The  simplest  assumption  which  will  approximate  the 
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complex  variation  of  the  surface  tension  law  with  bubble  growth  is  that 

S(R,c  ) is  a linear  function  of  R.  Then  S * (R  - R0)e,/R{>(n  - 1)  for 

nR  . If  one  puts  r * R/R  . then 
o o 


R_  ^ R £ nR,  and  S = tr  when  R 
o o 

the  surface  tension  law  can  be  written  in  the  form 


r - 1 


: n 


S = < 


(3) 


, n 


It  must  he  borne  in  mind  tha,t  the  surface  tension  law  defined  above  is 
in  a sr-nr-e  air  attempt  to  account  for  an  avenge  behavior  of  a large  number 
of  nuclei  of  many  possible  initial  sizes.  Compared  with  the  atomic  or 
molecular  scale,  the  nuclei  are  macroscopic  structures  and  hence  the  laws 
of  surface  tension  for  macroscopic  systems  are  applicable.  It  has  been 
argued  that  in  order  that  such  systems  may  exist  in  a stable  state  it  is 
necessary  to  add  a solid  phase  to  the  gas-water  system,  and  this  has  re- 
sulted in  the  introduction  of  another  parameter,  namely,  the  slope,  l/n, 
of  the  surface  tension  law.  To  ascribe  an  average  behavior  to  such  a large 
number  of  nuclei  with  many  possible  initial  sizes  is  reasonable  because 
experiments  show  that  the  zone  of  cavitation  numbers  at  which  cavitation 
starts  is  very  narrow.  Thus  it  is  possible  to  assign  a definite  value  of  the 
cavitation  number  for  the  inception  of  cavitation  for  a given  flow  situation. 


Primary  Parameters  and  Pres  sure -Time  Relationships 

The  problem  then  is  to  study  the  behavior  of  the  "average”  nucleus 
in  the  transient  pressure  regions  caused  by  the  flow  of  the  nuclei-containing 
water  around  submerged  bodies  of  various  sizes  at  various  flow  velocities. 
Thus,  the  relationships  between  several  parameters  must  be  studied  in 
order  to  obtain  approximate  quantitative  results  which  may  then  be  used 
to  guide  further  experimental  work.  In  particular,  one  must,  study  the  ef- 
fects upon  incipient  cavitation  of  bubble  air  content,  p , of  the  nucleus 
size,  Rq,  and  of  changes  in  the  slope  of  the  surface  tension  law  (l/n). 
Furthermore,  estimates  for  the  initial  radius  will  be  found  to  depend  upon 
the  nucleus  air  content  and  the  slope  of  the  surface  tension  law,  so  that 
all  of  these  factors  are  interrelated.  However,  it  will  be  possible  to 
make  some  theoretical  conclusions  about  the  effects  of  air  content  and 
surface  tension  law  slope.  To  test  such  findings  more  detailed  experimental 


results  are  required.  Unfortunately,  there  are  at  present  almost  no  re* 
liable  experimental  data  for  the  effect1  of  dissolved  air  content  U£on  the 
inception  of  cavitation.  Studies  are  tfow  under  way  for  determination  of 
the  effect  of  dissolved  air  content  on  incipient  cavitation. 

With  regard  to  pressure -time  relationships,  no  exact  account  will 
be  taken  of  the  many  pressure  distributions  which  can  arise  from  all  of 
the  various  body  shapes  that  one  may  wish  to  consider  in  a liquid  flow. 
However,  the  general  pressure  change,  first  falling  and  then  rising, 
common  to  all  such  bodies,  will  be  app-rr-^mated  by  two  parabola -like 
curves  joined  at  the  point  of  minimum  pressure.  The  term  "parabola- 
like" is  U3ed  because  the  actual  functions  used  will  be  parabolic  in  time 
and  the  resulting  spatial  pressure  distributions  will  differ  slightly  from 
the  true  parabolic  shape.  Two  curves  joined  at  the  minimum  pressure 
point  are  used  so  that  actual  pressure  distributions,  which  are  seldom 
symmetrical  about  the  minimum  pressure  point,  can  be  more  closely 
matched.  It  is  the  general  qualitative  behavior  of  the  many  possible 
pressure  distributions  that  is  of  importance  here.  No  attempt  will  be 
made  to  find  an  exact  forcing  function  for  any  specific  body  shape. 

Next,  certain  conditions  will  be  determined  which  will  enable  the 
"parabolic"  pressure  distribution  to  be  replaced  by  a dynamically  equiva- 
lent step  function  pressure  distribution.  The  detailed  numerical  integra- 
tions of  the  equation  of  motion  for  the  bubble  will  be  made  by  using  this 
xurtner  simplification.  Two  distinct  advantages  ate  gained  by  the  use  of 
the  step  function.  First,  it  turns  cut  that  for  a given  n,  or  surface  ten- 
sion law  slope,  if  the  parabolic  approximation  is  used,  a three  parameter 
family  of  solutions  for  Eq.  (1)  must  be  found.  However,  when  the  step 
function  pressure  distribution  is  used,  only  a two  parameter  family  of 
solutions  is  required.  Second,  the  use  of  the  equivalent  step  function 
allows  one  to  find  certain  important  relationships  between  the  primary 
parameters  involved  in  the  problem.  These  relationships  are  r.~t 
readily  determined  if  the  pressure  is  not  constant  during  some  time  inter- 
val. 

It  was  mentioned  above  that  the  only  connection  between  the  bubble 
growth  and  the  liquid  flow  over  the  body  is  the  choice  of  P(t)  in  the  equa- 
tion of  motion  for  the  bubble  wall.  For  the  step  function  pressure  distri- 
bution it  is  evident  that  P(t)  can  be  written  as 
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p and  V are  the  free  stream  static  pressure  and  velocity,  re  - 
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where 

spectively,  p is  the  liquid  density,  and  a is  the  pressure  coefficient  in 

the  low  pressure  region.  Fig.  1 
shows  the  step  function  pressure 
JL| Time,  t Uw. 

if  }'  is  reo'.*»***b»  rzd  tha, 
the  bubble  is  assumed  to  grow 

isotherm&lly,  and  if  account  is 

taken  of  Eqs.  (3)  and  (4),  the 

equation  of  motion  (Eq.  1 ) 

Fig.  1 becomes 
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p t > 0.  (5) 


The  initial  conditions  are  R(0)  = Rq  and  R(0)  = 0. 

Estimates  for  Initial  Effective  Bubble  Radius 

From  inspection  of  Eqs.  (3)  and  (5)  one  sees  that  the  forces  tending 

to  retard  the  bubble  growth  are  a maximum  at  R = nRQ.  If  the  bubble 

grows  in  such  a manner  that  R just  reaches  nRQ  with  aero  velocity,  it 

turns  out  that  the  time  required  for  such  growth  will  be  infinite  and 

cavitation  will  not  occur.  On  the  other  hand,  if  the  forces  reach  equi- 

librium  when  R = nR  he  motion  will  have  considerable  momentum,  and 

o 

the  bubble  growth  will  be  little  influenced  by  the  retarding  forces.  Cavi- 
tation will  then  be  well  started.  One  can  use  this  condition  of  force  equi- 
librium together  with  experimental  values  for  cavitating  flow  to  obtain  an 
estimate  of  the  initial  radius  Rq.  Thus,  setting  the  right  hand  side  of 
Eq.  (5)  equal  to  aero,  one  finds  that  for  R = nRQ 

. j 2 Pa  ^ 
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Fig.  2 ' ^>ates  for  the  initial  effective  bubble  radius  under 
various  dissolved  air  and  surface  tension  law  conditions 


where  K is  the  cavitation  number,  (pQ  - p^/(l/Z  pV^).  Experimental  value; 
of  K and  Vo  were  taken  from  the  data  for  incipient  cavitation  number  for 
the  Joukowski  hydrofoils/  and  curves  of  Rq  versus  n for  various  air  con- 
tents, p , were  calculated.  The  results  are  shown  in  Fig.  2.  (V  = 30  fps  = 
914  cm/sec,  K = 0.30,  or  = 70  dynes/cm,  a = 0.53).  Inspection  of  the 
Rq  vs.  n curves  shows  that  except  for  the  case  of  no  dissolved  air  in  the 
water,  the  variations  in  the  initial  radius  with  the  surface  tension  law  slope 
are  not  large.  Further,  since  Rq  is  given  by 


n(a  - K)(l/2  pVQ  ) + pa/n 
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an  underestimate  of  the  pressure  difference  (a  - K)  results  in  an  over- 
estimate of  Rq.  In  using  experimental  values  of  (a  - K)  for  incipient 
cavitatioit,  one  must  recognize  that  such  an  underestimate  of  (a  - K)  is 
being  made.  The  chief  value  of  the  calculations  is  that  they  show  that  an 
initial  radius  of  the  order  of  10  *cm  is  reasonable  (Fig.  2). 


Critical  Conditions  for  Cavitation 

As  was  mentioned  above,  the  forces  tending  to  retard  bubble  growth 
reach  a maximum  at  R = n R^.  However,  if  the  bubble  motion  acquires 
enough  momentum  before  reaching  this  point,  the  bubble  will  grov  •btougb 
this  critical  region  even  if  the  net  force  acting  on  the  bubble  in  this  region 
tends  to  collapse  it.  On  the  other  hand,  for  certain  values  of  the  coef- 
ficients in  the  equation  of  motion  (5),  the  small  radius  n Rq  will  not  be 
greatly  exceeded,  and  by  definition  the  condition  for  incipient  cavitation 
will  not  be  fulfilled.  In  this  section  relationships  will  be  derived  between 
these  coefficients  which  will  define  the  threshold  of  visible  cavitation. 


Before  deriving  the  relationships  between  the  coefficients,  the 
equation  of  motion  will  be  rewritten  in  dimensionless  form.  If  one  puts 

r = R(t)/RQ  and  r = t V 2<r  /p  Rq3  Eq.  (5)  becomes 
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where  the  "bubble  driving  parameter"  o is  given  by 
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and  the  "air  content  parameter"  Y is 
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The  initial  conditions  are  then 


r ( 0 ) = 1 , r o 

dx 
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where  the  initial  conditions,  Eq.  (9),  have  been  used  and  the  two  integrals 
have  been  matched  at  r * n.  Equation  (10)  simply  states  that  the  change 
in  the  kinetic  energy  during  the  expansion  is  equal  to  the  work  done  in  the 
course  of  the  motion.  If  one  lets  W(r)  represent  the  terms  on  the  right- 
hand  side  of  Eq.  10,  then  for  small  enough  values  of  the  bubble  driving 

parameter,  a ,•  the  function  W(r)  is 
very  near  a cubic  polynomial,  as 
shown  schematically  in  Fig.  3.  The 
value  of  r,  for  which  W has  its  iso- 
lated minimum,  will  be  called  r^. 

It  can  be  shown  that  for  n greater 
than  1 the  minimum  of  W occurs  at 
a value  of  greater  than  n.  Since 
the  kinetic  energy  is  zero  for  r = 1, 
it  will  again  be  zero  at  the  minimum 
of  W if  yfxnixi  ~ 0*  It  can  also  be 
shown  that  the  time  required  for  the 

bubble  to  grow  to  the  value  r^,  which  corresponds  to  a zero  minimum  of 
W,  is  logarithmically  infinite.  Thus,  if  the  bubble  driving  parameter,  o , 
and  the  bubble  air  content  parameter,  Y,  are  chosen  in  such  a way  that 
Wmin  = 0*  then  cavitation  will  not  occur  under  any  circumstance.  One  may 
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call  such  values  of  the  bubble  driving  parameter,  a , and  the  bubble  air 

content  parameter,  Y,  critical  values.  These  critical  values  are  denoted 

by  ac  and  Y , respectively.  A relationship  between  the  critical  values 

a and  Y will  now  be  found  from  the  conditions  W(r.)  * 0 and 
c c '1' 

dW(r^)/dr  a 0. 

From  Eq.(lO)  the  minimum  conditions  give 

o , ■■  ' • 1 r?  n2+  n + 1 

W(r,)  . -£(r,  - 1)  ♦ V*'l  ' ~ * — l 

and 

dW(ri)  2 Y 

- a r,  4 — - r.  = 0 . 

dr  rj 

Because  the  minimum  occurs  for  r^  > n,  only  the  part  of  the  inte- 
gral (10)  for  rkn  U required.  If  r^  is  regarded  as  a fixed  parameter, 
ac  and  Yc  are  given  in  the  parametric  form. 
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The  minimum  value  of  r.  corresponds  to  Y =0  and  the  root  of  0 a 

2 2 1 ^ 
r^  - ^r^  of  physical  significance  is  r^  = */Qc"  Substitution  of  this  re- 
sult for  r i into  ac  (Eq.  11)  gives 
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Equation  (12)  has  one  real  root  for  ac<  When  ac  is  real,  it  is  expedient 
to  solve  for  n in  terms  of  qc<  Since  n 1,  Eq.  (12)  gives  for  the  posi- 
tive root 


" = 2 


r 3 2 

8o  -3oS4 
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Fig.  4 - The  critical  conditions  for  the  occurrence  of  cavitation 

Taking  values  of  oc  from  0.1  to  1.0,  one  finds  the  corresponding  values 
of  n and  r^.  These  results  were  plotted.  For  various  fixed  values  of  the 
surface  tension  slope  number,  n,  the  values  of  r.  and  a for  Y = 0 

ICC 

were  taken  from  this  curve  and  then  larger  values  of  r ^ were  selected  to 
compute  and  Yc  from  Eqs.  (11).  The  results  of  these  computations  are 
graphically  presented  in  Fig.  4.  These  curves  of  the  critical  driving 
parameter,  a^,  as  a function  of  the  critical  air-content  parameter,  Yc  * 
are  the  loci  of  points  for  which  the  time  required  for  a bubble  to  grow 
from  its  initial  sine  to  a size  near  nRo*  (or  larger)  is  infinite.  For 
example,  if  for  a given  slope  (l/n)  of  the  surface  tension  law,  a value  of 


^Actually,  the  growth  is  from  r * 1 to  r=r^^n  with  r j 10n. 
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a greater  than  ac  is  chosen  at  a fixed  value  of  the  air  content  parameter, 

Y,  a finite  time  will  be  required  for  a bubble  to  reach  a definite  radius. 

. Further,  if  the  slope  of  the  surface  tension  curve  is  increased  (n  decreased) 

the  bubble  growth  time  increases  without  limit  as  a approaches  o.c> 

Numerical  Calculations  wi»h  the  Step  Function  Pressure  Law 

In  order  that  definite  numerical  calculations  can  be  mads,  values  of 

n and  Rq  must  be  chosen.  If  the  original  bubble  is  stabilized  on  only  a 

portion  of  an  unwetted  solid  particle,  it  seems  plausible  that  the  bubble 

must  grow  to  several  times  its  initial  size  ^rfore  the  c -jnjsion  law 

reaches  its  full  value  <r  . For  the  present  calctua;.UM»  a value  of  5 will 

be  taken  for  the  parameter,  n,  in  the  surface  tension  law.  If  one  takes 

account  of  the  overestimate  contained  in  the  curves  of  R versus  n 

° -4 

(Fig.  2),  a value  for  the  initial  effective  radius,  Rq,  of  2 x 10  cm. 
seems  reasonable. 

In  accordance  with  the  definition  of  incipient  cavitation  given  above, 

the  bubble  must  grow  to  a macroscopic  size  of  R = 1 millimeter.  The 

total  range  of  r,  (s  R(t)/R  ) will  then  be  from  1 to  500.  The  problem  of 
• ^ 

determining  the  condition  for  incipient  cavitation  is  now  reduced  to  the 
following  question.  What  time  interval,  as  a function  of  the  bubble- 
driving parameter  a and  the  bubble  air-content  parameter  Y is  required 
for  the  bubble  to  reach  a radius  which  is  500  times  its  initial  value  ? To 
answer  this  question  one  must  find  a two-parameter  family  of  solutions 
of  £q.  (10)  of  the  form  t = F(a,Y).  The  dimensionless  time  parameter,  t, 
is  a function  of  the  body  length,  free  stream  velocity,  and  cavitation  number. 

The  air  content  parameter,  Y,  will  be  given  and  a , the  bubble -driving 
parameter,  is  a function  of  the  cavitation  number  and  free  stream  velocity. 

Then  if  t = F(  a , Y)  is  known,  a trial  and  error  procedure  will  give  the  in- 
cipient cavitation  number  when  the  free  stream  velocity  is  known. 

Equation  (10)  is  easily  reduced  to  quadratures,  so  that  the  dimension- 
less time  T is  given  by 


“16" 


The  two  integrals  arise  because  cf  the  change  in  the  surface  tension  law 
at  r = n = 5.  Since  the  variations  in  t are  required  for  small  changes  in 
a and  Y,  approximate  representations  for  Eq.  (13)  will  be  limited  in  value 
for  those  cases  where  a is  near  the  critical  bubble  “driving  parameter  qc 
In  any  event,  numerical  integrations  have  to  be  precise  for  the  first  inte- 
gral in  Eq.  (13). 

For  the  range  of  the  air  content  parameter,  V ,from  0 to  1.8,  with 
corresponding  values  of  the  bubble -driving  parameter,  a,  from  0. 145  to 
0.230,  ninety-eight  integrations  of  Eq.  (13)  were  carried  out  on  an  I.  B.M. 
calculator.  From  these  results  the  functional  relationship,  t = F(o,  V), 
represented  by  Eq.  (13)  was  plotted  as  in  Fig.  5.  Tabulated  data,  from 
which  the  figure  was  made,  are  given  in  Table  I. 


After  the  above  trial  and  error  calculations  for  incipient  cavitation 
number  versus  velocity  have  been  made  for  different  values  of  the  air 
content  parameter  Y,  comparisons  between  experimental  results  and  the 
calculations  for  the  given  shape  could  be  used  to  obtain  a measure  of  the 
success  with  which  the  surface  tension  law  parameter,  n,  was  chosen. 

It  will  be  noticed  from  Fig.  4,  that  the  slopes  of  the  curves  for  critical 
bubble -driving  parameter,  ac,  versus  the  critical  bubble  air  content 
parameter,  Yc,  are  steepest  for  low  values  ol  the  surface  tension  law 
parameter.  Thus,  the  cavitation  number  for  incipient  cavitation  will  be 
smaller  for  steeper  slopes  of  the  surface  tension  law,  S(R,  a),  than  it  will 
for  a surface  tension  law  with  more  gradual  slope  (large  n).  If  sufficiently 
reliable  data  are  available,  a comparison  with  the  calculated  curves  for 
incipient  cavitation  will  suggest  a better  choice  for  the  effective  surface 
tension  law  parameter,  n.  Because  there  are  almost  no  reliable  data  at 
this  time,  these  alterations  of  the  calculations  will  not  be  made  here. 
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Propertles  of  the  "Parabolic"  Pressure  Function 


Aa  noted  previously,  the  step  function  results  will  be  related  to  the 

growth  of  a bubble  under  a more  realistic  pressure -time  forcing  function. 

The  important  general  features  of  many  such  actual  forcing  iunctions  can 

be  qualitatively  represented  by  two  functions  which  are  parabolic  in  time. 

One  function  which  decreases  from  C = 0 to  C » -b  and  another  func- 

P P 

tion  which  increases  from  C * -b  to  C * 0 will  be  employed.  Two 

P P 

functions  are  used  because  the  pressure  distributions  around  actual  bodies 
are  seldom  syrw  ,-trical  in  the  streamwise  direction.  Fig.  6 shows  the 
parabolic  pressure  distribution.  The  required  results  for  the  "parabolic" 
pressure  distribution  will  now  be  obtained. 


Decreasing  part,  0«S&  0^^  -b: 

The  equation  of  a time  parabola  which  decreases  from  =•  0 to 
-fc  in  a time  interval  T • may  be  written  as 


(14) 


where  t is  any  time  in  the  interval  0 t ^ T 
(Fig.  6).  From  the  Bernoulli  equation  dz/dt  = 
V = VQ  Vl  - Cp(x)  , where  a is  the  distance 
along  the  stream  line  adjacent  to  the  body. 
Under  the  conditions  x = x(t),  with  dx/dt  \ 0 
one  may  write  dx/dt  = VQ  Vl  * Cp(t).  If 
x{0)  = 0 and  x(T)  = ip  one  obtains 


Fig.  6 - The  time -parabolic 


or 


(15) 
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When  Cp  = - K,  that  is,  the  pressure  coefficient  on  the  body  is  equal  to 

minus  the  cavitation  number  fo?  the  flow,  the  pressure  in  the  water  around 

the  bubble  is  equal  to  the  vapor  pressure  of  the  water.  Inspection  of  Eq. 

(5)  shows  that  the  bubble  can  not  start  to  expand  until  it  has  reached  a point 

on  the  model  where  C = - K.  The  region  for  which  C ^ - K will  be 

P P 

called  the  region  of  "favorable  environment"  for  bubble  growth.  Thus  the 
time  that  is  of  concern  here  is  the  time  during  which  the  bubble  is  in  a 
favorable  environment.  The  value  of  (t/T)  for  which  Cp  = - K is  ob- 
tained directly  from  Eq.  (14).  If  tj  is  the  required  value  of  the  time, 
then 


b-K 


T 


but  the  time  interval  Atj  spent  by  the  bubble  in  the  favorable  environ- 
ment is  (T  — tj).  Hence  AtjaT^(b  — K)/b,  or  substituting  for  T one 
finds 


2 yfb  - K 


(y/b  + (1  f b)  sin"1  yJjjK 


The  relationship  x = x(t)  can  be  found  by  suitably  altering  the  limits  of 
integration  in  the  integral  leading  to  Eq.  (15).  If  Eq.  (15)  is  solved  for  lj 
and  divided  into  the  result  for  x = x(t),  one  obtains 


ti./l  + b n t.  , 1 + b . -1,,  t. 
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Increasing  part,  - b ^ Cp  ^ C : 

For  the  portion  of  the  pressure  distribution  downstream  from  the 
minimum  pressure  point,  one  can  write 


CpW 


b (t_ 

(m  - 1 VT 


n2  2 1 
- 1 ) - (m  - 1) 
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where  Thu  a,  Cp(T)  * - b and  Cp(mT)  * 0.  In  a manner 

similar  to  that  used  in  obtaining  Eq.  (15),  the  analogous  result 


l !, 


(m  - 1)  T = 


{VF*h*-jV£)* 


(19) 


is  obtained.  Here  is  the  distance  along  the  stream  line  next  to  the 
body  from  the  point  where  Cp  = - b to  the  point  where  again  Cp  = 0.  and 
T is  given  by  Eq.  (15).  The  time  t^  at  which  Cp  = - K is  obtained  from 
£q.  (18) 


— - 1 - (m-1) 
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b-K 


The  period  spent  by  the  bubble  in  the  favorable  environment  for  this  por- 
tion of  the  pressure  distribution  is  given  by  a t^  - (t^  - T).  Substitution 
for  (m  - 1)T  from  Eq.  (19)  yeilds 


At, 
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As  before,  the  space -time  relationship  can  be  found  in  a manner  similar 
to  that  employed  in  deriving  Eq.  (17).  Thus, 
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where  x is  measured  along  the  stream  line  from  the  point  for  which 
Cp  * -b.  From  Eq.  (15)  the  expression  for  T can  be  put. into  Eq.  (19)  to 
find  m.  This  substitution  gives 


(m  - 1)  = -i 

*1 


(22) 
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The  "parabolic"  curves  are  compared  with  the  Joukowski  and  the  "Hemis- 
phere" experimental  pressure  distributions  in  Figs.  8 and  10.  Equations 
(16)  and  (20)  can  be  added  to  find  the  total  time,  at  = At^  + ht^,  that  a 
bubble  spends  in  the  whole  region  of  favorable  environment,  and  one  finds 

2(1. +i  ) yV-  K 

i At  - i — t 

It  is  t'nis  value  of  the  time  At  which  determines  the  dimensionless  time  t. 
For  a given  pressure  distribution,  1^,  and  b are  known  quantities  and 
Vo#  the  free  stream  velocity,  will  be  given  so  that  K is  the  only  unknown 
factor  required  to  find  At  or  t. 

Dynamic  Equivalence  between  the  Step  and  "Parabolic"  Pressure  Laws 

The  preceding  preliminary  calculations  now  make  it  possible  to  ap- 
proximate the  relationship  which  must  obtain  between  the  idealized  step 
function  pressure  distribution  and  the  more  realistic  "parabola -like" 
pressure  distribution  if  the  essential  features  of  the  bubble  growths  are 
to  be  the  same  in  both  cases.  The  matching  will  consist  of  two  parts. 

First,  the  time  spent  by  the  bubble  in  the  low  pressure  region  will  be  taken 
to  be  the  same  for  both  the  step  function  and  the  parabolic  pressure  laws. 
Second,  the  equation  of  motion  (1)  will  be  integrated  in  closed  form  by  ap- 
proximate means  for  both  the  step  function  and  the  parabolic  pressure 
laws.  The  two  approximate  bubble  histories  will  be  said  to  be  dynamically 
equivalent  when  the  total  bubble  growths  under  the  two  pressure  laws  are 
1 equal.  These  two  conditions  will  result  in  approximations  of  the  required 

relationships  between  the  step  and  parabolic  pressure  functions. 

The  first  point  of  comparison  for  the  two  pressure  laws  is  the  re- 
lationship between  the  time  intervals  during  which  the  nucleus  is  exposed 
to  the  favorable  environment  in  each  case.  It  will  be  specified  that  the 
time  spent  by  the  bubble  in  the  favorable  invironment  shall  be  the  same 
* for  both  the  step  function  pressure  law  and  the  "parabolic"  pressure  dis- 

tribution. In  both  instances  the  free  stream  velocity,  VQ,  will  be  identi- 
cal. From  2q.  (4),  the  pressure  coefficient  for  the  pressure  step  is 


1 


22 


for  t <C  0 , cr  x < 0 . 
for  t 0,  or  x ^ 0, 


(24) 


W 

If  the  bubble  is  assumed  to  grow  in  an  interval  0 £x  ^ X,  then  the  same 
reasoning  which  gave  £q.  (15)  gives  t * — — — ~ — , where  At  is  the 


total  time  during  which  the  bubble  is  exposed  to  the  low  pressure  and  X is 
the  distance  along  a streamline  traversed  by  the  bubble  in  the  time  At. 
Equating  this  result  to  the  A t of  Eq.  (23)  gives 


X 2(i.  + I,)  Vb“  K 

= * — j=r  (25) 

Vl+a  VT  ♦(Ub).ia-VT^ 

If  a relationship  between  a and  b can  be  found,  then  X will  be  given  in 
terms  of  known  quantities. 

The  second  point  of  comparison  is  the  requirement  that  the  total 
bubble  growth  in  the  time  interval  of  favorable  environment  must  be  the 
same  for  the  two  pressure  functions.  This  requirement,  which  will  result 
in  a relationship  between  a and  b,  will  be  called  dynamical  equivalence. 

The  a priori  derivation  of  an  exact  condition  for  the  dynamical  equivalence 
of  the  two  cases  is  formidable  in  view  of  the  fact  that  to  obtain  such  a pre- 
cise result  the  equation  of  motion  (1)  must  be  integrated  in  exact  analytical 

form.  It  will  be  worthwhile  to  use  an  approximate  method  of  integration 

8 

due  to  M.  S.  Plesset.  Plesset's  method,  while  lacking  rigorous  justi- 
fication, has  been  shown  to  give  results  in  close  agreement  with  precise 
numerical  results  if  the  right-hand  side  of  the  dynamical  Eq.  (1)  is  a 
function  which  increases  exponentially  with  time.  In  this  instance, 

Plesset's  method  approximated  the  exact  numerical  results  within  0.3% 
for  sufficiently  large  values  of  the  bubble  growth  time.  On  the  other  hand, 
if  the  pressure  is  constant,  it  can  be  shown  that  this  approximation  will 
give  the  final  bubble  radius  with  an  error  of  about  5%. * Thus,  if  the  right- 
hand  side  of  the  equation  of  motion,  Eq.  (1),  is  a monotonic  function  of 

*This  comparison  approximates  the  effect  of  surface  tension  for  bubbles 
growing  through  the  specified  size  range  because,  when  surface  tension  and 
air  content  are  negligible,  Plesset's  method  is  exact.  This  comparison  was 
made  for  bubbles  with  surface  tension  and  initial  internal  air  pressure. 
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timc,  it  is  evident  that  a suitable  approximation  for  the  final  bubble  size 
will  result.  Plesset's  method  is  based  on  the  supposition  that  the  right- 
hand  side  of  the  equation  of  motion  is  a monotonic  function  of  time.  It  is 
unfortunate  that  this  condition  will  not  be  entirely  satisfied  here.  After 
obtaining  the  required  approximate  result,  it  will  be  used  to  define  a pair 
of  dynamically  equivalent  pressure  distributions  so  that  numerical  inte- 
grations can  be  performed  to  test  the  reasonableness  of  the  approximate 
method. 

The  equations  of  motion  applied  to  each  case  differ  only  in  the  form 
of  the  forcing  function  on  the  right-hand  side,  because  in  one  case  a step 
function  is  taken  while  in  the  other  case  a "parabolic"  dependence  is  as- 
sumed. Thus,  it  seems  reasonable  that  in  comparing  the  behavior  in  the 
two  cases,  one  can  discard  the  surface  tension  and  air  content  terms. 

This  approximation  may  be  partially  justified  by  recognition  of  the  fact 
that  both  the  surface  tension  and  the  air  content  terms  decrease  very 
rapidly  as  the  bubble  grows.  Thus  Eq.  (1)  can  be  written  in  the  form 

RR  + 1 R2  = g(t)  , (26) 

2 


where  g(t)  is  given  for  the  step  function  case  by 


g(0 


V0Z(-K)  . 

VQ2  (a  - K)  , 


t < 0 , 
t 0 ; 


and  for  the  parabolic  case  by 


g(t) 


b-K 

b 


with 

o c o - 4)  * 


(27) 


(28) 


Equation  (28)  accounts  for  only  the  portion  of  the  parabolic  pressure 
curve  lying  upstream  from  the  minimum  pressure  point  which  is  in  the  re- 
gion of  favorable  environment.  The  growth  in  the  second  portion  of  the 
parabolic  pressure  curve  will  be  accounted  for  by  approximating  the  slope 
of  the  growth  curve  at  the  minimum  pressure  point  and  multiplying  this 


} 


slope  by  the  time  interval  of  favorable  environment  (Eq.  20). 

Pleaset's  approximate  integration  method  it  baaed  upon  the  condi- 
tion '.hat  in  Eq.  (26)  g(t)  ia  a function  that  doe  a not  permit  the  acceleration 
of  the  bubble  wall  ft  to  become  negative;  then 


R ft  = g(t)  - | ft*^0. 


Consequently: 


kZ4Ll  g(t)  . or 


7 S(t) 


Integration  of  (30)  gives 


R - R 


jf  • 


where  R ia  the  radiua  at  the  initial  time  t . Now  define  a function  $(t) 
o o 

such  that 

R - R0  = 9(t)^^/|g(e)  d£  =<p(t)I(t)  . (3i 

Because  of  Eq.  (31),  the  function  (p(t)  is  of  bounded  variation: 

N 

o s£  <p(t)  ^ i . (3: 

Equation  (32)  is  a formal  solution  to  Eq.  (26).  Accordingly,  substitution 
of  (32)  into  Eq.  (26)  gives  the  exact  differential  equation, 

(R  + *1)  ($1+291  + 9!)  + |(<pl  + i<p)2  = |i2  • (3' 

O b b 

Liquation  (34)  will  now  be  approximated  with  the  assumptions  that 

9 I + 24>  I « 9 i , 

9 1 <5C  9 I , 


and  R <££  <p  I . 

o 

A first  approximation  for  <p  is  then  given  algebraically  by  Eq.  ;34). 


-25  “ 


Solving  for  9 yields, 


<P  S3 


(35) 


This  result,  together  with  Eq.  (32),  gives  the  required  approximate  inte- 
gration. In  order  that  the  parabolic  case  can  be  treated  by  means  of  the 
linear  extension  outlined  above,  the  derivative  dP.(t)/dt  must  be  found.  If  one 
takes  account  of  the  inequalities  preceding  Eq.  (35)  the  derivative  is 
given  by  R sts  <*  t.  From  the  results  following  F.q.  (19)  dt^  = (t^  - T)  = 

(m  - 1)  T V(b  - K)/b  so  that  the  linear  growth  to  be  added  is 


R(T)  Atz 


<T  i (m  - 1)  T 


Hence  the  total  growth  in  the  parabolic  case  is 


R - R <P  (T) 
o ' 


I(T)  ♦ I(T)(m  - 1)  T 


b - K 


(36) 


Substituting  the  equation  for  g(t)  (Eq.  28)  into  the  integral  for  i(t), 
one  finds, 

t 

'T 


i(t) . VT-^vo2b[^ii-(1-u)23<,u ; 


V3 


b-K 


or  putting  y = (1  - u) 


I(t)  = - TV, 


' V 3b 


‘0 


-1 )J±- 

T V b-K 


V 1 — y^  dy  = - 


2 


J(b  -K); 
V 3b 


+ sin  * (1 


"T + (l  "T*  (i  “ 
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He nee, 


i (t)  = jl  t v a / fc  r.ftf 
4 ° V 3b 


From 


i(t)  = v. 


fb  -K 


b-K\  T/ 


t \2  ~ V 

] and  I (t)  - — — 


K) 


t\  , / b-K 


T V1  “t) 


it  follows  that, 


i(T)  =V0\/^i^  and  I(T)  = 0 . 


From  Eq.  (38)  and  Eq.  (35),  <p  (T)  *s  1 so  that  R - Rq,  Eq.  (36), 


becomes 


R-R  JbTV  li+m-1 
o o \ . 


(b.-isr 


for  the  parabolic  case. 


i(t>  = 


f \ wo  (a  - K>  dt’ 


Because  the  time  spent  in  the  favorable  environment  is  the  same  in  this 
case  as  it  is  in  the  parabolic  case,  the  time  interval  for  Eq.  (40)  is  equal 
to  the  At  of  Eq.  (23).  It  is  easy  to  see  that  At  = mT  yj( b - K)/b  so 


I ( At)  = m VQ  T y . . (41) 

Inspection  of  Eq.  (40)  shows  that  Y(t)  = 0 so  that  <p  (t)  - 1.  For  this  case 
Eq.  (32)  alone  suffices  to  give  the  approximate  integration  required,  and 
hence 


R - R a,  m V T A >=■ tlllZ tlL  (42) 

° V 3 b 

The  requirement  for  dynamical  equivalence  can  be  applied  to  Eqs.  (39)  and 
(42)  to  give 
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m \/a  - K **  (i+  m - 1)  Ajb  - K . (43) 

Thus,  the  relative  shapes  of  the  step  and  the  parabolic  pressure  distribu- 
tions are  now  specified  by 


2(4l  + l2) 

Vb  + (1  4-o)  sin  1 


(25) 


/ tt/4  + m - 1\  

V a - K *s  f < JVb  “ K , (43) 

and 

(m  - 1)  = f2Aj  (22) 

The  above  matching  conditions  were  checked  for  one  case  of  incipi- 
ent cavitation  on  the  two  inch  Joukowski  hydrofoil.  From  the  trial  and 
error  calculations  presented  below  (p.  29  ).  it  was  found  that  when  the  free 
stream  velocity  is  16  meters  per  second  the  incipient  cavitation  number, 

K,  is  0.403,  and  hence  the  bubble -driving  parameter,  o,  is  0.190  and 
the  dimensionless  time.T,  is  1400.  The  equation  of  motion  (1)  was  inte- 
grated for  the  case  where  the  air  content  parameter,  V,  is  1.0  and  the 
P(t)  term  is  composed  of  the  parabolic  laws  (16)  and  (20).  The  numerical 
solution  gives  a bubble  diameter  of  2/3  rr«tr».  at  the  end  of  the  specified 
time  interval.  If  the  matching  had  been  exact  the  final  bubble  diameter 
would  have  been  2 mm.  The  two  results  are  comparable,  well  within  an 
order  of  magnitude.  Hence,  one  may  cor. elude  that  the  matching  condition 
(Eq.  43)  between  the  step  and  parabolic  pressure  laws  will  give  reasonable 
results  for  the  present  application. 

Scaling  Laws  for  Incipient  Cavitation  at  High  Free  Stream  Velocity  (Vn— » oo) 

It  is  now  possible  to  obtain  a relationship  for  cavitation  scale  effect 

for  very  large  values  of  the  velocity,  V . As  the  free  stream  velocity, 

V , becomes  very  large,  the  bubble -driving  parameter,  q , (see  Eq.  7) 

° 2 

Increases  as  V From  Eq.  (23)  it  is  clear  that  the  time,  At  (and  hence 
o t — — 

T since  x = At  y 2o"  /p  Rq  ) decreases  as  l/V  . Hence  for  a » 1, 


*28* 


Eq.  (13)  can  be  approximated  by 


T «« 


m 


where  a = (a  -K) 


RoPvfc 

o o 

4e’ 


r/4  + m-l\  R pV 


m 


(b-K, 


This  result  for  a follows  from  Eqs.  (7)  and  (43).  From  Eq.  (23)  the  rela* 
tion, 


T = 


\jb  - K 


vo  (V*  + <» +b> sin'1  Vtte-  ) 


2<r 

• 

PR' 


will  be  used  for  t in  Eq.  (44).  For  geometrically  similar  bodies,  ( Jj  + 1^) 
will  be  some  constant  multiplied  by  the  characteristic  body  length  L.  Using 
this  fact  and  substituting  for  t and  a into  Eq.  (44),  one  finds 


L V b-K 


V(b“K> 


(45) 


where  C is  a constant.  Suppose  now  that  a given  body  running  at  very  high 
velocity  has  a characteristic  length  Lq  and  is  found  to  have  an  incipient 
cavitation  number  K0-  Thrr  if  another  r.imilar  body  in  high  speed  flow  has 
length  L and  incipient  cavitation  number  K,  the  relationship  (b  - KQ)/(b  - K) 
s l/Lq  must  hold.  This  last  result  can  be  written  in  the  form 


subject  to  the  condition  that  V , the  free  stream  velocity,  is  very  high. 
Just  how  high  the  velocity  Vq  must  be  is  determined  from  more  detailed 
calculations  below.  A plot  of  K/Kq  versus  l/Lq  is  given  in  Fig.  7 for 
values  of  b/KQ  from  1.0  to  2.0. 
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Fig.  7 - The  scaling  condition  for  incipient  cavitation  for 
very  large  free  stream  velocities  (Vo— 1 ► oo) 

Calculati  r for  Incipient  Cavitation  on  the  Joukowski  Hydrofoils 

Calculations  for  incipient  cavitation  on  the  Joukowski  hydrofoils  will 
now  be  made  by  means  of  the  procedures  developed  above.  The  relation- 
ships which  must  be  used  for  the  numerical  work  are 

v2  n _ ,,  2 

( f + m - i 

a * 


/J  + m-lV  R PV  * 

- (b-K)  ; 

\ m ) 


A* 


(m  - 1)  = l2/fj 

*0*!  + l2)  ^b-K 


T * 


(/7 


♦ (1  ♦ b)  sin 


'V?k> 


omsoaiH  tmuMUNt  otmoctt  alom  boot. 


and  a graph  of  the  function  t * F(o,  y), 
Fig.  5.  In  order  to  find  m,  the  experi- 
mental pressure  distributions  for  the 
three  Joukowsky  hydrofoils  were  cor- 
rected for  tunnel  blockage  (1),  and  a 
single  curve  was  faired  through  the 
points,  as  shown  in  Fig.  Q.  It  was  ob- 
served experimentally  that  the  first 
visible  cavitation  on  the  hydrofoils 
occurred  at  about  20%  of  the  chord,  c, 
from  the  point  where  the  pressure 


L. 


-30- 


distribution  Cp  is  first  zero.  In  tn  effort  to  account  for  this  experimental 
fact,  the  distance  i ^ was  taken  to  be  equal  to  0.25  c - Ij.  The  length 
was  determined  by  measuring  the  distance  from  the  minimum  pressure 
point  to  that  point  where  the  pressure  coefficient,  Cp,  is  zero  nearest  to 
the  leading  edge  of  the  hydrofoil.  These  measurements  gave  ^/l. 

14/11  so  that  (m  - 1)  * 1.  27.  These  geometric  relationships  and  the  ap- 
proximating parabolic  pressure  distribution  are  shown  in  Fig.  8.  If  one 

•A  3 

takes  «r  = 70  dynes/cm,  Rq  = 2 x 10*'  cm,  pa  1 gram/cm  , b = 0.53, 
the  equation  for  the  dimensionless  time  is 

T = 1.28  x l<^  y]0.S3  - K , 

vo 

and  for  the  bubble -driving  parameter, 

c = 0.581  x 10‘6  (0.53  - K)  VQ2  . 

By  means  of.  the  graph  for  t = F(o,  Y)»  with  Y = 1.0  , trial  and  error 
calculations  were  made  to  find  K for  values  of  the  free  stream  velocity 
VQ  of  from  900  cm/sec  to  3000  cm/sec  for  each  of  the  three  chord  lengths,' 
2 in.  (5.08  cm),  4 in.  and  8 in.  The  calculations  were  extended  to  infinite 
velocity  by  means  of  Eq.  (45)  (for  this  case  C a 550  cm).  For  the  case 
c = 4 in. , the  calculations  were  repeated  for  zero  air  content  (y  = 0).  From 
these  calculations,  curves  of  incipient  cavitation  number  K versus  the  ratio 
c/VQ  are  given  in  Fig.  9.  The  calculated  results  are  also  compared  with 
experiment  in  Fig-.  9<  The  experimental  scatter  is  too  great  to  draw  any 
definite  conclusions  as  to  the  applicability  of  the  theoretical  results.  The 
calculated  curves  show  the  tendency  for  the  cavitation  number  to  increase 
with  velocity  for  all  scales,  and,  except  for  the  higher  velocities,  no  scale 
effects  are  found  at  constant  velocity.  The  data  do  not  conclusively  prove 
or  disprove  these  theoretical  findings. 


Y = 1 corresponds  to  an  initial  bubble  air  pressure  of  700  millibars 
or  a dissolved  air  concentration  of  70%,  the  saturation  concentration  for 
an  air  pressure  of  1 atmosphere  above  the  water. 


Calculations  for  Incipient  Cavitation  on  the  "Hemispheres" 


Calculations  were  made  for  the  inception  of  cavitation  on  right  cylin- 
drical bodies  with  hemispherical  noses.  The  same  body  sizes  were  se- 
lected as  were  used,  by  R r.srmeen  in  his  experiments  (see  Ref.  1). 
Kermeen  found  that  the  dissolved  air  concentration  in  the  water  during  his 
experiments  was  about  one -half  of  the  saturation  value  for  an  air  pressure 
of  one  atmosphere.  The  corresponding  value  of  the  air  content  parameter, 
y,  was  taken  to  be  0.  7. 

a 

The  experiments  of  Rou^e  show  that  b, 
the  absolute  value  of  the  minimum  pressure 
coefficient  on  the  body,  increases  from  about 
0.  65  at  a Reynolds  number  of  0.  6 x 10®  to 
0. 74  when  the  Reynolds  number  is  greater 

5 

than  or  equal  to  1. 2 x 10  . This  change  in 
the  minimum  pressure,  b,  was  accounted  for 
in  the  determinations  of  the  dimensionless 
time,  Ti  (Eq.  23)  and  of  the  bubble -driving 
parameter,  o , (Eqs.  7 and  43).  The  "para- 
bolic"  pressure  distribution  was  fitted  to  the 
experimental  pressure  distribution  by  selecting  for  the  pressure  distribu- 
tion asymmetry  factor,  m,  a value  of  1.45.  Figure  1C  compares  the  ex- 
perimental and  the  "parabolic"  pressure  distributions  when  the  Reynolds 
number  is  in  the  supercritical  regime  (Re  ^ l.Zx  10®).  Trial  and  error 
calculations  were  then  performed  in  the  same  fashion  as  for  the  Joukowski 
hydrofoils.  The  calculations  are  compared  with  Kermeen's  data  in  Fig.  11. 
It  will  be  seen  that  the  calculated  values  of  incipient  cavitation  number  are 
in  poor  agreement  with  experiment.  The  disparity  between  the  experi- 
mental and  the  calculated  results  increases. as  the  sixe  of  the  body  de- 
creases. There  are  two  possible  reasons  for  the  disagreement.  Either 
the  mathematical  approximations  used  in  developing  the  theory  are  too 
great,  or  some  physically  important  effect  (such  as  that  due  to  the  bound- 
ary layer)  has  been  overlooked.  It  may  well  be  that  both  of  these  factors 
contribute  to  the  large  differences  between  theory  and  experiment. 

It  will  be  noticed  in  Fig.  11  that  for  the  2,  1 -l/8  and  l/4  inch  models 
the  experimental  data  for  incipient  cavitation  number  seem  to  show  little 


Fig.  10 
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or  no  dependence  upon  the  free  stream  velocity,  VQ  , when  VQ  is  greater 
than  or  eaual  to  80  ft  per  sec.  This  apparent  independence  of  the  incipi- 
ent cavitation  number  with  free  stream  velocity  corresponds  to  the  theo- 
retical conditions  leading  to  Eq.  (46).  Equation  (46)  relates  the  incipient 
cavitation  number  to  the  characteristic  body  length  of  similar  bodies 

unde:  the  condition  that  the  free 
stream  velocity  is  very  high.  The 
incipient  cavitation  number  under 
this  condition  is  said  to  have  at- 
tained its  limiting  value.  If  it  is 
supposed  that  the  experimental 
data,  for  the  highest  free  stream 
velocities  measured,  approximate 
such  limiting  values  of  the  incipi- 
ent cavitation  number,  it  is  reason- 
able to  compare  the  experimental 
points  for  this  case  with  Eq.  (46). 
The  data  for  incipient  cavitation  r^mber  on  each  model  at  the  highest  test 
velocity  were  av-raged  and  plotted  against  model  diameter  in  Fig.  12. 
Equation  (46)  wac  .!,Cu  fitted  to  the  experimental  points  twice  at  diameters 
of  l/4  in.  and  2 in.  , respectively.  The  two  curves  are  shown  in  Fig.  12. 

It  will  be  seen  that  the  agreement  between  experiment  and  theory  is  best 
when  the  calcinations  are  fitted  to  me  smallest  model.  The  free  stream 
velocity  in  all  ca  .es  was  at  least  80  't  > . • Perhaps  data  for  even 

higher  values  of  Vc  would  give  better  agreement  when  the  theory  is  com- 
pared with  experiment.  More  experiments  must  be  made  to  test  this  point. 
The  theoretical  trends  (Figs  9 0 r.d  II)  jeem  to  indicate  that  the  limiting 
value  for  the  incipient  cavitation  number  will  be  attained  for  higher  veloci- 
ties, VQ  , for  larje  bodies  than  for  small  models.  If  this  trend  is  correct 
it  means  that  most  water  tunnels  cannot  reach  high  enough  speeds  to  make 
use  of  this  theoretical  approximation  when  models  with  low  pressure  re- 
gions larger  than  that  of  the  two-inch  diameter  hemisphere  are  tested. 

It  will  be  recalled  that  in  deriving  Eq.  (46),  as  with  all  other  theo- 
retical results,  two  basic  approximations  were  made.  First,  the  pressure 
was  estimated  from  experimental  results  and  Bernoulli's  Law.  Second, 
the  time  during  which  the  nucleus  is  exposed  to  a low  pressure  condition 


Fig.  ii 


1 


-35- 


where  bubble  growth  can  occur  was  estimated  from  the  body  length  and  the 
liquid  velocity  along  the  body  just  outside  the  boundary  layer.  The  appli- 
cation of  Bernoulli's  Law  should  give  quite  accurate  results.  One  must 
conclude,  therefore,  that  the  time  estimates  are  not  good  enough  for  the 
present  purposes.  Furthermore,  if  one  traces  through  the  calculation 
procedure  for  the  inception  of  cavitation,  he  will  find  that  underestimating 
the  time  results  in  an  underestimate  of  the  cavitation  number  for  incipient 
cavitation.  The  present  calculations  do  not  account  for  the  possibility  that 
at  least  a portion  of  the  bubble  growth  occurs  in  the  boundary  layer  where 
the  time  available  for  growth  wouiri  be  greater  than  the  ti  ne  estimates 
used  here.  Thus,  before  one  can  make  predictions  for  the  inception  of 
cavitation  with  confidence,  the  boundary  layer's  effect  '-n  growing  cavita- 
tion bubbles  must  be  understood. 
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PART  III 
COKCi  TJSIONS* 

It  is  found  by  experiment  that  < avitation  on  submerged  bodies  in  a 
rectilinear  flow  shows  decided  changes  in  its  development  as  the  free 
stream  velocity  and  body  size  are  changed.  It  is  found  that  under  certain 
conditions  the  incipient  cavitation  number  depends  on  both  body  scale  and 
free  stream  velocity.  Of  these  two  effects,  that  due  to  changes  in  free 
stream  velocity  is  the  greater,  and  most  of  the  monotonic  increase  in 
incipient  cavitation  number  with  free  stream  velocity  occurs  below  a free 
stream  velocity  of  60  fps.  In  the  velocity  range  between  30  and  60  fps  the 
cavitation  number  changes  by  about  25%  of  its  maximum  value.  For  ve- 
locities above  60  fps  the  incipient  cavitation  number  increases  by  less  than 
10%  of  its  maximum  value.  At  constant  velocity  the  incipient  cavitation 
number  increases  with  scale.  However,  more  work  must  be  done  before 
truly  quantitative  results  can  be  given  for  the  scale  effect  at  constant 
velocity. 

It  was  found  that  one  cannot  obtain  useful  scaling  laws  for  limited 
cavitation  from  similarity  arguments.  In  fact,  the  similarity  calculations 
show  that  one  may  expect  limited  cavitation  to  be  affected  by  both  body  size 
and  free  stream  velocity.  One  is  then  led  to  conclude  that  the  cavitation 
number,  K,  is  not  the  only  significant  parameter  required  to  define  the 
limited  cavitation  flow  regime.  Both  experiment  and  analysis  indicate 
that  one  must  specify  the  model  size,  free  stream  velocity,  dissolved  air 
content  and  the  cavitation  number  if  he  wishes  to  describe  completely  an 
experimental  situation  for  an  immersed  body  of  specified  shape.  The 
condition  of  the  body  surface  should  also  be  specified,  but  this  effect  has 
not  been  '-onsidered  in  this  paper.  Although  it  is  convenient  to  use  the 
cavitation  number,  K,  for  theoretical  calculations,  one  would  be  closer 
to  the  physical  relationships  if,  instead  of  the  cavitation  parameter,  he 
specified  the  liquid  vapor  pressure  and  the  free  stream  static  pressure. 

From  the  lack  of  agreement  between  experiment  and  the  greatly 
simplified  theoretical  calculations  presented  in  this  report,  one  can 

^These  conclusions  summarize  the  results  of  both  the  preliminary 
veport(I)  and  the  v/ork  presented  in  this  paper. 
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conclude  that  the  neglect  of  the  role  of  the  bovindary  layer  is  not  justified 
for  precise  predictions  of  the  effects  of  model  scale  and  free  stream  ve- 
locity upon  the  inception  of  cavitation.  Further,  the  use  of  an  arbitrary 
linal  radius  instead  of  a specified  maximum  radius  for  the  required  babble 
growth  may  have  introduced  an  extraneous  variability  in  the  calculations. 

On  the  other  hand,  the  qualitative  agreemeni  between  the  experimental  and 
calculated  trends  for  incipient  cavitation  number  versus  free  stream  ve- 
locity and  body  size,  substantiates  the  basic  premise  that  scale  effects  in 
limited  cavitation  arise  because  of  transient  pressure  effects  on  the  nuclei 
in  the  liquid  flow. 

A combination  of  simple  theoretical  considerations  and  some  of  the 

expo rin'icntal  data  for  incipient  cavitation  gives  estimates  for  the  effective 

initial  nucleus  radius.  It  is  found  that  for  the  water  used  in  the  experiments, 

-4  -3 

the  calculated  effective  initial  radius  is  between  10  and  10  centimeters. 

For  the  limiting  case  where  the  free  stream  velocity  approaches 
infinity,  it  is  found  that  the  cavitation  number,  K,  for  incipient  cavitation 
does  not  approach  b,  the  absolute  value  of  the  minimum  pressure  coef- 
ficient on  the  submerged  bedy.  Instead,  the  limiting  value  of  the  incipient 
cavitation  number  approaches  a value  less  than  b,  and  the  difference, 
b — K,  is  found  to  be  inversely  proportional  to  the  characteristic  body 
length,  L.  The  theory  closely  approximates  the  experimental  results  if 
the  equation  is  fitted  to  the  data  for  the  smallest  model.  More  experi- 
mental work  is  needed  to  demonstrate  the  complete  validity  of  this  result. 
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